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Abstract. We present new polarisation observations at 13 cm, acquired using the Westerbork Synthesis Radio 
Telescope (WSRT), of 65 sources, from the B3-VLA sample of Compact Steep-Spectrum sources. These new 
data are combined with our VLA polarisation data, at 3.6, 6 and, 21 cm, presented in a previous paper. Due 
to the multi-channel frequency capabilities of the WSRT, these new 13 cm observations enable a more reliable 
determination of integrated Rotation Measures, and of depolarisation behaviour with wavelength. The new data 
are inconsistent with the depolarisation models that we used earlier, and we propose an alternative model which 
seems to work properly. We also revise our previous model for the external Faraday screen, and its dependence 
on the source redshift. 
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1. Introduction 

Compact Steep-Spectrum (CSS) sources and GHz-Peaked 
Spectrum (GPS) sources, because of their small size, are 
fully embedded in the Interstellar Medium (ISM) of the 
host galaxy. Their radio properties are affected by the 
properties of the ISM, unlike large-size radio sources which 
extend well beyond the optical dimensions of their host 
galaxies, and can therefore be used to probe the condi- 
tions of the ISM. Polarisation characteristics are a useful 
tool for this purpose. 

The measurable quantities from polarisation observa- 
tions are: 

— The Faraday Rotation Measure (RM), defined as the 
slope of a linear fit of the polarisation angles, as a 
function of A 2 : 

X (\) = RM \ 2 + X (0). 

If the medium is homogeneous (both its density and 
magnetic field) , or if the inhomogeneities are resolved 
by the observing beam, xW is strictly proportional 
to A 2 at all wavelengths, and RM = k f Bundi (this 
quantity is also called Faraday Depth), where Bt\ is 
the component of the magnetic field along the line of 
sight, n is the electron density of the medium, I is 
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the geometrical depth of the medium along the line of 
sight (los), and k is a constant. If the medium is un- 
resolved, or partially resolved by the observing beam, 
RM changes from point to point across the source, dif- 
ferent contributions of polarised radiation are rotated 
differently, and x(A) can deviate, more or less strongly, 
from the A 2 -linear law. 
— The change of the fractional polarisation m\ with A, 
described by the quantity DP(X) = m\/m . If the 
medium is uniform, or if its inhomogeneities are well 
resolved by the observing beam, DP(X) = 1. In an 
inhomogeneous medium, the variations of RM from 
point to point will produce changes in m\ versus A. 
DP(X) will exhibit different behaviours with A, de- 
pending on the properties of the inhomogeneities in- 
side the medium. 

Several models h ave been developed to inter pret the 
Faraday effe cts (e.g. Burn 1966t iTribble I fl991r . see also 
iLaing lfl984l for an excellent and concise review). The be- 
haviour of the polarisation angle, x(A), and fractional po- 
larisation, m\, as a function of A 2 , enables the "average" 
Faraday Rotation Measure, and, from the "screen mod- 
els", its dispersion, <trm, to be determined. Using these 
data, we can obtain information on the density distribu- 
tion of the ISM that surrounds the radio source, on its 
dumpiness, and on both the ordered and tangled compo- 
nents of the magnetic field. 
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Polarisation studies of CSS and GPS source sam- 
ples have been conduct e d by several authors (e.g . 



van Brc ugel et al. 



Staiighcll ini et al. 



1984 lAkuior fc Garrington I [l995; 



19981 : IPeck fc Tavlorl l2000h . They 



have found that GPS sources are almost unpolarised, while 
the larger-size CSS sources can show large RMs and/or 

large depolarisations, as a f unction of A 2 . 

In a previous paper (jFanti et al. 1 12004 , hereafter 



Paper IV), we used "low resolution" polarisa tion measures 
at 8.5 and 4.9 GHz 3.6 and 6 cm) from iFanti et al. I 
(|200ll hereaf ter Paper I), and at 1 .4 GHz 21 cm) from 



the NVSS (jCondon et al. I Il998f) . to derive the RM 



and the depolarisation properties of a comp lete sample 
of CSSs, the B3-VLA CSS sample (|Paper | . Our main 
results were: 



1. 



In general, the total source depolarisation was found 
to follow either the Burn or the Tribble model. 
In > 50% of the cases, the integrated \W follows the 
A 2 -linear law from 3.6 to 20 cm. The derived RMs have 
values of up to a few hundred radm~ 2 . After subtrac- 
tion of the Galactic Rotation, and correction for the 
source redshift, z, we found that « 20% of the sources 
have intrinsic RMs of up to 1000 radm -2 . 
There is a wavelength-dependent characteristic scale 
(from sa 3 kpc at 3.6 cm to > 6 kpc at 21 cm), below 
whi ch radio sources are a lmost totally depolarised (see 
also lCotton et al. |[2003h . 

crm increases with redshift; a similar, but less signifi- 
cant, dependenc e was suggested for RM (Figs. 13 and 
14 in IPaper IVl ). 



To explain our results, we proposed a simple model, 
based on Faraday effects, with an appropriate spatial dis- 
tribution of the a mbient gas density and magnetic field. 

The results of IPaper IVl . however, were based on po- 
larisation data with a non-optimal wavelength coverage, 
because the gap in A 2 , between 6 and 20 cm, is too wide. 
In a number of cases there were remaining ambiguities in 
both RM and <7rm- In fact, for the few sources for which a 
polarisation measurement was available at the intermedi- 
ate wavelength of 11 cm dKlein et al. 1120031 ) the initial 
model was not suppo rted in a number of cases (see Figs. 4 
and 10 in IPaper Pvt ). 

To improve our polarisation information we per- 
formed new polarisation observations at 13 cm, using the 
Westerbork Synthesis Radio Telescope (WSRT) for 65 ra- 
dio sources, 58 of which were detected in polarisation at 
one or more of the three available VLA frequencies. The 
remaining 7 unpolarised objects were observed as control 
sources. 

The new observations fill the large gap in A 2 between 
the 6 and 21 cm VLA data. They allow us to improve 
the reliability of the RMs significantly by reducing the 
ambiguities, and constrain the depolarisation behaviour 
as a function of A 2 . 

Section [5] provides a short description of the B3-VLA 
CSS sample and of the previous polarisation data (3.6, 6, 



and 21 cm, VLA), and presents the selection criteria for 
the WSRT sub-sample. 

Sections [3] describes the new WSRT polarisation ob- 
servations at 13 cm, the data reduction strategy, and the 
derived results. 

Section [3] summarizes the polarisation status of the 
WSRT sub-sample. 

Section [5] presents the results on Rotation Measure 
(RM) and its dispersion (ctrm)- 

Section [6] revisits the model of the ambient magneto- 
ionic medium. 

Section [7] provides our conclusions. 

Appendix [Al contains the data table and comments on 
individual sources. 

Appendix [B] describes a simple two-polarised- 
component model, which has been applied to a minority 
of radio sources. 



2. The WSRT sample 

The sources discussed i n this paper were sele cted from the 
B3-VLA CSS sample (|Vigotti et al. Ill989h . described in 
Paper J . which consists of 87 CSSs/GPSs (three of which 
do not have polarisation data) with flux density > 0.8 Jy 
at 408 MHz, with projected Linear Sizes (L5]0in the range 
0.2 k- 1 < LS (kpc) < 20 h' 1 . Their radio luminosity, at 
the selection frequency, is Po.4GHz > 10 26 hr 2 WHz -1 . 
The sources were observed using the VLA in A configu- 
ration, a t 6 and 3.6 cm in total intensity and polarisa- 
tion (see IPaper J ). A detailed desc ription of the polarisa- 
tion data reduction was provided in lPaper" IV . In addition, 
polarisation data at 2 1 cm are available from the NVSS 
1998h . 



(Condon et al. 



A sub-sample of 65 sources, hereafter referred to as 
the "WSRT sub-sample" , was observed using the WSRT 
at 13 cm. At this wavelength, most sources, according to 
the previous VLA total polarisation measurements, were 
expected to have a polarised flux density Sp J> 1 mJy. 
Seven sources, undetected in polarisation at 3.6, 6 and 
21 cm, were observed as a control sample. 

The 19 B3-VLA CSS sources that were not observed 
at 13 cm, were either unpolarised, or strongly depolarised 
at 6 cm and 20 cm. 

The WSRT sub-sample includes 45 out of the 54 B3- 
VLA CSS sources of Linear Sizes larger than 2.5 kpc. Of 
the missing sources, 2 are unpolarised at all frequencies, 
6 are polarised at 3.6 cm only, and one, polarised at 3.6 
and 6 cm, does not meet the 1-mJy selection criterion. 



1 In this paper we have kept Ho = 100 h kms -1 Mpc -1 , 
and qo = 0.5 for consistency with previous papers. We 
have also used the Concordance Cosmology, with Ho = 
71 h kms" 1 Mpc -1 , Q M = 0.27, fi V ac = 0.73, and found that 
the results discussed in Sect. [6] remain largely unchanged. 
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3. The Polarisation Data 

3.1. The WSRT Observations 

The observations were carried out in November 2004, 
using the WSRT, in dual circular polarisation mode, at 
the mean frequency of 2263 MHz (w 13 cm). Eight 
Intermediate Frequencies (IFs) were used, each one di- 
vided into 64 identical frequency channels. 

Using the full available bandwidth of 128 MHz, the 
expected root mean square (rms) noise in the I, U, Q 
Stokes parameters (Mi,u,q)i for a 20-minute integration, 
is ~ 0.14 mJy/beam. The 8 x 64 = 512 independent fre- 
quency channels enable radio interferences to be removed 
reliably, which helps to bring the noise level close to its 
theoretical value. 

Thanks to its sensitivity and high resolution (~ 8" at 
13 cm), the WSRT is at present the only instrument, in 
the northern hemisphere, that can detect such low-levels 
of polarised flux-density, at 13 cm. The Effelsberg tele- 
scope, for example, is confusion-limited in polarisation at 
a level of ~ 0.5 mJy ( Klein et al. Il20"03h at the closeby 
wavelength of 11 cm (beam size sa 4.'3). 

We observed each target source in snapshot mode at 
three well-spaced hour-angles of approximately (0 h , ±4 h ), 
for a total integration time of ~ 20 minutes. Flux-density 
and phase-calibrator sources were observed, on average, 
every 4 hours, for 15 minutes. 



3.2. Data reduction 

All data reduction (editing, calibration, imaging and 
analysis) was performed with the NRAO package AIPS 
(Astronomical Image Processing System). Interference 
spikes were removed using the task UVLIN. 

3.2.1. Calibration 

3C 286 was used as a primary calibrator, for flux-density, 
phase, bandpass and polarisation angle. Secondary cali- 
brators were 3C 147 and CTD 093 (known to be unpo- 
larised at 13 cm). The flux-density calibration uncertainty 
was ~ 1%. The flux density scale is within 3% of that of 
Baars et all (|l978l ). 



To determine the residual instrumental polarisation 
(D-term), the task LPCAL was applied to the observa- 
tion of CTD 093. The measured value was approximately 
0.1% of the source flux density (Si), which is consistent 
with the results for the 7 unpolarised sources of the sam- 
ple (Sect. [5]). After the D-term calibration, an arbitrary 
offset in the polarisation angle x remains, which was de- 
termined by using integrated measurements of the source 
3C286. Each IF was corrected separately. We were un- 
able to obtain a good calibration for IF1 and IF7. The 
data reduction was, therefore, based on 6 IFs out of 8 for 
polarisation data, while total intensity data were derived 
using all 8 IFs. 



Based on the r.m.s. of the polarisation angles obtained 
for each IF and each scan of 3C 286, we estimate that the 
polarisation-angle calibration is accurate to within 1?4. 

3.2.2. Imaging 

The three snapshots of each source were combined (task 
IMAGR), to produce two-dimensional "dirty images" of 
1024 x 1024 pixels (~ 43' x 43') for the Stokes parameters I, 
Q, U, and occasionally V. This size was usually sufficiently 
wide to identify and remove all confusing field sources oc- 
casionally present in polarisation. We then cleaned the 
images down to the theoretical noise level. 

At the WSRT resolution (« 8"), all sources were un- 
resolved. We derived total-band flux densities (S 1 , S Q , 
S v ) for the Stokes parameters I, Q, U by fitting a bidi- 
mensional Gaussian to the brightness distribution (task 
IMFIT), and setting the search boxes for Q and U about 
the source, as visible on the total intensity image (I). 
When the Q and/or U signal were too weak for reliable val- 
ues to be produced using IMFIT, we used instead IMEAN, 
which integrates the surface brightness inside the box. 
Using these measurements, we derived the polarised flux 
density Sp - 



S l 



S 2 ,. 



A noise-dependent statistical correction was applied 
to the polarised flux-densities to c o rrect for polarisation 
bias dWardle fc Kronberg 1974 ; ISimmons fc Stewart I 



1985HCondon et al. II 19981). The de-biased polarised flux- 



density is S£ b = ^/Sj 

of Sp (see Sect. l3.2.'31 . According to IWardle fc Kronberg 



jip, where /xp is the noise error 



(|1974| ). this formula is appropriate for Sp/^ip ^1-2. The 
fractional polarisation, m — Sp b /Si, was also obtained. 

Thirty-three sources were detected in polarisation at 
13 cm, at levels > 3ixp, using the total bandwidth. 

3.2.3. Noise estimate and errors 

The pixel histograms in "empty" regions of the I, Q 
and U images, and the pixel statistics of the V im- 
ages, are approximately Gaussians, and are in agreement 
with each other. We can therefore assume that the rms 
is a reliable noise estimate. The typical r.m.s. error is 
w 0.15 mJy/beam for all Stokes parameters, in agreement 
with expectations. 

We adopted noise errors appropriate to each individual 
source, and we used the statistics of the residuals provided 
by AIPS, after the Gaussian fit, to derive the values of 
fj,i, /xq, /Ltij- We assume that this procedure takes account 
of the fit quality, and the possible confusion by residual 
sidelobes from field sources. The distribution of the indi- 
vidual noise errors is approximately Gaussian, ranging in 
value between 0.08 mJy/beam and 0.24 mJy/beam, with 
a peak at « 0.15 mJy/beam, in agreement with the above 
estimate. 

For /iq and fj,y, when IMEAN had to be used, the 
noise errors were computed using the above r.m.s. ~ 
0.15 mJy/beam, scaled by the square-root of the num- 
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ber of beam areas included in the search area. The total 
noise error, /^p, was computed by taking into account the 
different values of /xq and The uncertainties of the 
flux-density calibration, and of the residual instrumental 
polarisation, were quadratically added to the noise error, 
to derive the final errors of Si and Sp h (a\ and crp). The 
error of m was computed using error propagation: 




The noise error of x was also computed, using error 
propagation: 



1 

2S\ 



2 V ( S QW) 2 + (SuMQ) 2 



According to IWardle fc Kronberg I (|l974f ). this formula 
is correct only for m > 2cr m . 

The total error, tr x , was computed by quadratically 
adding the calibration error of « 1?4 to the noise error. 

3.2.4. In-band Polarisation Data 

For 24 sources with signal-to-noise ratios J> fO, we mea- 
sured the polarisation parameters for each IF (in-band 
data). For all but two sources, we detected signal in the 
individual IFs, with a signal-to- noise ratio J> 5. 

For the determination of the flux densities of the Stokes 
parameters in individual IFs, we used the same procedure 
as for the total band (Sect. I3.2.2[) . The errors of Sq and 
S\j were estimated from the r.m.s. of the six Sp of the in- 
dividual IFs. The average value was crp = 0.5 mJy/beam, 
rather than the expected 0.4 mJy/beam. This value is 
consistent with the r.m.s. of the in-band polarisation an- 
gles, computed for sources that show no significant gradi- 
ent of the angle across the band. We therefore empirically 
adopted the value crp = 0.5 mJy/beam as the error of the 
polarised flux density Sp of each IF. 

No use was made of the in-band polarised flux density, 
while the x from the individual IFs were linearly fitted, 
and the derived in-band RM13 used as a guide, to solve 
the mr ambiguities of the polarisation angles, when deter- 
mining the total RM. 

3.2.5. "Recalibration" of the Polarisation position 
angle 

The polarisation position angle was computed initially as 
X = \ tan _1 (5u/S , Q). When we compared these angles 
with those obtained by interp olation at 13 cm, with a 
A 2 law, of the VLA data from iPaper IVI . we found large 
disagreements. At this point we decided to use, as inter- 
nal calibrators, the 8 sources B3 0110+401, B3 0213+412, 
B3 0754+396, B3 0800+472, B3 0805+406, B3 0955+390, 
B3 1220+408, arid B3 1 343+386, which were also observed 
by iKlein et al.l (|2003l ), using the Effelsberg telescope, at 
the close wavelength A = 11 cm. The polarisation angles 
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Fig. 1. Comparison between Xwsrt and XEf f ci S bcr g - 
Errorbars of WSRT data do not usually exceed the plotted 
symbol size. The sources in the plot are: 3C 286 (primary 
calibrator; large cross at [33,33]) and 8 sources obse rved 
using the Effelsberg telescope by IKlein et al.l (|2003l ) (see 
text). The r.m.s. about the A 2 -linear fit is « 10°, compat- 
ible with the errors. 



measured with the WSRT, Xwsrt, and those measured at 
Effelsberg, XEffci S bc Ig , are plotted in Fig. [T] 

The two sets of angles are clearly related to each 
other by the relation XEffcisbc-g = —Xwsrt + 66°. Obviously 
XEffdsbcrg and Xwsrt coincide for 3C286, because the po- 
larisation angles were calibrated using this source. 

Given the AIPS definition for crossed circular polari- 
sation, RL = Q +i\J, LR = Q -iU, i.e. U = i(LR - RL), 
Q = RL + LR, (i — \J — 1), our result implies that there 
is a swap between RL and LR. 

We corrected empirically all of the WSRT angles, using 
the above relation. In Fig. [21 we show a few examples 
of sources before and after the swap of the cross-hand 
visibilities on the 13 cm angles. We also plot the in-band 
angles (see Sect. 13.2^4"!) . and confirm the correctness of our 
approach. 

4. The WSRT Complete Polarisation Data Set 

The polarisation status of the 65 sources, observed using 
the WSRT at 13 cm, is provided below: 

— 26 sources have been detected at > 3 op, at all four 
wavelengths (labelled PI in Table |A~T|) : 

— 4 sources have been detected at > 3 crp, at 3.6, 6, and 



13 cm (labelled P2 in Table [OJ; 

6 sources have been detected at > 3 crp, 

6 cm (labelled P3 in Table [OJ; 



at 3.6 and 
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Fig. 2. Examples of sources with \i3 before (o) and after (•) the empirical correction of Sect. 13.2.51 The original 
in-band 13-cm angles are largely displaced from the general x(A) behaviour (x ) and also generally have a wrong slope. 
Moreover, the trends before the empirical correction have a A 2 slope with opposite sign to that after the correction. 
The corrected ones have the slope in agreement with the general RM. The last panel (B3 1055+404 A) shows how the 
13-cm data helped in solving the nir ambiguity of the fit obtained with three wavelengths only (□ represent a l-ir 
ambiguity in \)- 



5 sources have been detected at > 3 er P , at 3.6 cm only polarisation is <^ 10% for RM ^ 400 radm 2 , while at 



(labelled P4 in Table EI}; 

— 8 sources have been detected > 3 op, at two or three 
non-contiguous wavelengths (labelled P5 in Table lA~T} ; 
only 3 of these have been detected at 13 cm; 

— 9 sources have been detected at 21 cm only (labelled 
P6 in Table EI}; 

— the 7 remaining sources are undetected in polarisation 
at all wavelengths and are labelled NP (not polarised) 
in Table EH 

We remark that bandwidth depolarisation is not a prob- 
lem at 3.6 and 6 cm, unless very high Rotation Measures 
(> 10 4 radm -2 ) are present. At 13 cm, bandwidth de- 



21 cm it can be large (> 40% for RM > 200 radm -2 ; see 



Condon et al71ll998h 



The polarisation data are presented in Table IA.1I We 
provide: the total flux density (S13), the fractional polari- 
sation (TO13) with the corresponding error, the polarisation 
position angle (X13, defined within ±90°), the polarisa- 
tion parameters derived in Sect. i.e. Rotation Measure 
(RM) and Rotation Measure dispersion (orm), both ob- 
served (obs) and in the source frame (sf), and the intrinsic 
fractional polarisation (mo), and Covering Factor (here- 
after f c ) (see Sections 15.21 and 15. 3p . To provide all of 
the data used in this paper we include the redshift, ei- 



6 



Rossetti et al.: B3-VLA CSSs: polarisation 



ther spec troscopic or photometric from K or R band data 
(iPaper l|), and the source projected total Linear Size (taken 
from Paper ivh . When no redshift was available, z = 1.05 
was assumecQ We provide notes on individual sources, in 
Appendix [Al 

We emphasize that these data are integrated over the 
entire source. All sources are unresolved at both 13 cm 
(WSRT) and 21 cm (NVSS) wavelengths. In contrast, 
the majority of sources were resolved, or partially re- 
solved, by the VLA at 3.6 and 6 cm. Hence, at these short 
wavelengths, the polarisation parameters were vectorially 
added over the source extension. 



5. Results 

5.1. The "Cotton Effect" at 13 cm 



Cotton et al. ( 20031 ) found that, at 21 cm, the fractional 
polarisation of radio sources with LS < 6 kpc, is very 
small, typically < 0.4%. For LS « 6 kpc the fractional 
polarisation rises abruptly to a median value of approxi 



Paper IV 



mately 1%, reaching values of up to 4-5%. In 
we extended the analysis of fractional polarisation versus 
Linear Size to shorter wavelengths, and found the same 
effect at 3.6 cm and 6 cm. The critical size, defined vi- 
sually, that separates polarised from unpolarised sources, 
was estimated to be ~ 2.5 kpc at 3.6 cm, and « 4 kpc at 
6 cm. 

In Fig. we show a plot similar to those discussed 
above, for the data at 13 cm presented in this paper. 
The critical scale for the drop in fractional polarisation is 
about 5 kpc. 



10 



~ 6 - 
B 4 



10 
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15 



20 



Fig. 3. The "Cotton Effect" at 13 cm. We note that most 
of the 19 unpolarised sources (not observed by us; see 
Sect. [2} would fall at small LS. 



2 As discussed in IPaper fl , this is the average of the 
spectroscopically-determined redshifts for the objects of the 
B3-VLA sample, not detectable at the limit of the POSS plates, 
which were later identified using deeper observations. 



In Sect. 16.11 we interpret the abrupt change in the 
distribution of m\, as a function of LS, in terms of a de- 
pendence of <7RM,sf (see Sect. 15. 2p on LS. 

5.2. Modelling the Depolarisation of individual sources 

The analysis of the fractional polarisation, as a function 
of A 2 , indicates that the majority of sources have a frac- 
tional polarisation that decreases with increasing wave- 
length. This implies that an inhomogencous magneto-ionic 
medium, a Faraday curtain, is present in front of the 
source. The effects of such a curtain depend on the prop- 
erties of the inhomogeneities in the ma gnetized m edium, 
which are often referred to as "cells" . In IPaper IVI. we in - 
terp reted the data u sing the models of both Burn I (1966) 
and iTribblel (|l99llh In each model, the screen completely 
co vers the source (/ c = 1). 
" (Il966h 



Burn 



assumes that the "cells" are much 
smaller than the beam size, and produce a random dis- 
tribution of RMs, across the source, which are zero on 
average and have a dispersion crm- The fractional polar- 
isation follows a Gaussian law (in A 2 ): 



(1) 



where too is the intrinsic (A = 0) fractional polarisation. 
We performed a simple Monte Carlo analysis that shows 
that the model is applicable up to ctrmA 2 ~ 1 rad, pro- 
vided that the cell area is much smaller than a hundredth 
of the entire source area. It is well known, however, that 
Eq. Q] predicts a fractional polarisation that is too low 
at long- wavelengths, to be able to reproduce the obser- 
yations. This discrepancy is reduced in Tribble's model 
(|Tribble Ill99lh . Also this model assumes that the RM is 
randomly distributed, and has a dispersion ctrm, but in 
additi on it intro duces a broad distribution of "cell" sizes. 

In IPaper IV , we fitted the three- wavelength VLA po- 
larisation data with either model, solving for the long- 
wavelength fractional polarisations, which are often too 
high for the Burn model. We now have data for a fourth 
wavelength, between 6 cm and 21 cm, that enables far 
tighter model fits to be achieved (see Fig.0|). 

Using the new 13-cm data, the number of sources with 
polarisation data that can be fitted using the Gaussian 
Burn mo del, is reduced with resp ect to the findings in 
Paper IVI . We also find that fits in IPaper IVI . made with 
the Tribble model, are generally in bad agreement with 
13-cm data. As a matter of fact, for a large fraction of 
sources, the fractional polarisation drops quickly between 
3.6 cm and 6—13 cm, and remains approximately constant. 
This is inconsistent with the predictions of either model 
described above (see e.g. plot for the source B3 0754+396 
in Fig. [4|. In such cases, the following empirical modifica- 
tion of Eq. [T] 



m 



/ c e" 



+ (1 - fc) 



(2) 



is more successful in reproducing the data. 
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Fig. 4. Examples of fits from Paper ivl (3 wavelengths 
only) with Burn's (Eq. [TJ dashed line), Tribblc's (solid 
line) and a partial coverage (Eq. [2l dotted line) models 
which use also the 13 cm data. The top panel shows that 
with the old VLA data B3 0039+412 could have been fit- 
ted with a Burn model, while the addition of the 13 cm 
data indicates that Eq. Ogives a much better fit. Note that 
below 13 cm also Tribble's model could have been satis- 
factory within the measurement errors. In the lower panel 
the only acceptable fit for B3 075 4+396 is with Eq.[l The 

(l200l . 



Klein et al. 



For 



data point at A = 11 cm is from 
A < 11 cm Tribble's model would have been satisfactory 
as well. 



The obvious interpretation of Eq. [2] is that, if a source 
is only partially covered by the screen, a fraction (1 — f c ) 
of the source radiation emerges non-depolarised, and 
maintains a constant level of fractional polarisation at 
long wavelengths. We remark that ctrm only concerns the 
source "covered fraction" . At this stage, we however con- 
sider the introduction of f c to be a mathematical artifice 
that improves the fit. We discuss the possible origins of 
this "partial coverage" in Sect. 16.81 



We stress that, for 2cr^ M A 4 <C 1 rad 2 , the three mod- 
els provide similar results and differences appear only at 
longer wavelengths, ctrm derived in the short- wavelength 
regime is therefore largely model-independent. 



We applied Eq.[2]to the 36 sources in the WSRT sub- 
sample (labelled PI, P2, PS in Sect. [TJ) that were well 
detected at 3.6 cm and 6 cm, using the 13 cm and 21 cm 
fractional polarisations or upper limits as well. A large 
fraction of these sources show a flattening in fractional 
polarisation (as, e.g., in Fig. [4j) , which requires that f c < 
1. The remaining sources do not show the flattening of 
m\ within our wavelength range and can be fitted using 
Eq. [TJ However, for homogeneity, we also fitted their data 
using Eq.[2j In these cases, we provide a lower limit to the 
Covering Factor, and an upper limit to the corresponding 
ctrm- In < 50% of the radio sources, f c > 0.9; in a number 
of sources, however, f c can be as low as 0.3. Figure [5] 
shows that the fraction of sources with small f c is larger 
for LS ^ 6 kpc. 
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Fig. 5. Covering Factor f c vs LS. Open symbols represent 
optically-unidentified radio sources. Their LS have been 
computed by assuming that z = 1.05 (see note on page [6]). 

Three sources (B3 0110+401, B31025+390B, 
B3 1216+402) out of these 36, show an oscillatory 
behaviour in their fractional polarisation, even with 
repolarisation at long wavelengths. This behaviour could 
be due to the beating of sub-components with different 
Rotation Measures (sec Appendix [B]). The data for 
B3 0110+401 were modelled in this way, and a good fit 
was achieved for both m\ and x(A) (see Fig. [6]). For the 
other two sources, the fit is more poorly constrained (see 
notes to Table [A~Tj) . 

In addition to the above 36 sources, another 5 were de- 
tected only at 3.6 cm (P4). We applied the Gaussian Burn 
model (f c w 1) to these source data, and derived lower 
limits to crRM,obs, which are consistent with the sources 
being strongly depolarised at A <; 6 cm. 

Data for another 8 sources, detected at two or three 
non-contiguous frequencies (PS), were not well fitted by 
either Eq. [TJ or [2] Several of them might have oscillations 
as a function of A 2 in the fractional polarisation. For these 
data, we applied the model described in Appendix [Bl with 
some success, (see notes to Table IA.1[) . Because of the 
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Fig. 6. B3 0110+401. Left Panels: m\ and x(A) data. Note that X21 is largely off the A 2 -linear fit of the higher 
frequency data. Right Panels: Fit of m\ and x(A) with the two-component-model of Appendix [Bj The dotted curve 
in the modelled m\ takes into account a little depolarisation of the two components. The dashed straight line in the 
x(A) plot represents the average RM. 



limited amount of data, the parameters may however be 
weakly constrained. 

Of the 9 sources found to be polarised only at 21 cm 
(P6), one source (B3 2302+402) has a fractional polarisa- 
tion much larger than the upper limits at the other wave- 
lengths. The remaining 8 have upper limits at 3.6 cm, 
6 cm, and 13 cm close to mzi, and are therefore very 
little polarised at all frequencies, consistent with a large 
ctrm, and an f c slightly less than unity. 

Also the 7 sources that are unpolarised at all fre- 
quencies (NP), are likely to have a large o"RM,obs (j£ 
500 radrn" 2 ). 

Table I A. 1 1 provides the best-fit Rotation Measure dis- 
persions in the observer's frame (<JRM,obs) and in the 
source frame (cxRM.sf)- The latter are computed using 
the redshift, either spectroscopic or photometric, from 
Table EU Formal errors are typically < 20%. 

From the 45 sources of the WSRT sub-sample for which 
depolarisation parameters (or limits) could be compute, 



38 have 2.5 < LS(kpc) < 20. Their intrinsic degree of 
polarisation has a median value of approximately 4.3%, 
an r.m.s. of 2.5%, and a distribution tail exten ding up to 

"Tvl these 



Paper 



12% (see Fig. [7]). As already noted in 
intrinsic fractional polarisations do not differ significantly 
from those of radio sources of tens to hundreds kpc sizes. 

5.3. Rotation Measures 

We used the E-vector position angles at 13 cm fy 13 ) and 
those at 3.6 cm, 6 cm, and 21 cm, reported in lPaper IVl . 
to derive the Faraday Rotation Measure, RM, by weighted 
linear interpolation of the data, with the A 2 -linear law 
x(A) = xo+A 2 i?M. We are aware that the assumption of a 
A 2 -linear law could be incorrect. As discussed in Sect. 15.21 
the presence of a few polarised sub-components with dif- 
ferent RM would produce modulations of the x(A) versus 
A 2 -linear relation (see Appendix |B|) . A A 2 -linear fit would 
provide an average RM that corresponds approximately 
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Fig. 7. Distribution of the intrinsic polarisation, mo, for 
the WSRT sub-sample. 

to that of the component with the highest polarised flux- 
density. Individual data points may then deviate signifi- 
cantly from an average A 2 -linear law (e.g., Fig. [SI top-right 
panel). 

The fitting procedure was applied to the 32 sources 
that were detected in polarisation at 13 cm, with Sp > 
lOcrp, and have at least two additional detections at the 
3 up level. As stated in Sect. 13.2.41 the 24 sources with 



S P > 



10a p were analysed by means of the individual 



IFs. Twenty-two of these sources were detected at a level 
of more than 5erp in the individual IFs, and we were 
able to derive the in-band Rotation Measure, RM13, with 
typical uncertainties in the range of 25 - 70 radm -2 . 
Then, we used RM 13 as a guide to resolve possible wr- 
ambiguities at the other wavelengths. For four sources 
(B3 0034+444 in Fig. ® B3 0137+401, B3 0814+441, and 
B3 0930+389) the first fitted i?Mhad to be changed dras- 
tically, to reach agreement with RM13. For two addi- 
tional sources (B3 1216+402 and B3 1350+432) the dis- 
agreement between RM and RM13 remains; however this 
is probably due to substructures in polarisation that pro- 
duce a modulation of x(A), over the A 2 -linear law (see 
the individual notes for these sources). Another source 
(B3 0120+405) shows strong disagreement between RM 
and RM13. Although we have been unable to develop a 
model to explain this discrepancy, we suppose that the ex- 
planation proposed for the previous two sources may also 
apply to this one. The value of RM provided in Table lAT] 
in this case, was derived by fitting the individual IFs, to- 
gether with data at other wavelengths. 

The observed RMs (uncorrected for the Galactic 
Faraday Rotation, i?Af Q b s ), are provided in Table [A~T1 A 
code "a" , marks the sources for which a A 2 -linear fit pro- 
vides a good chi-square (probability > 0.05 of being ex- 
ceeded because of random fluctuations). Only 12 sources 
out of 32 are found to be in this class. The remaining 
sources have A 2 -linear fits of lower chi-square quality. 
In a number of cases, this can be attributed to a single 
discrepant point; the exclusion of this point would signifi- 
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Fig. 8. B3 0034+444. Example of how the in-band polar- 
isation angles help in choosing the RM. The positive RM 
(dashed line) is clearly ruled out. With the total band X13 
only, instead, both fits would have been acceptable. 

cantly improve the chi-squared fit, with moderate changes 
in RM. We suspect that, in a number of cases, we observe 
modulations, about the A 2 -linear law, caused by polarised 
sub-components of a source that have different RMs. For 
this reason, we excluded no data point from the fit. 

The formal errors of RM are small, typically 
10 radm -2 . The actual uncertainties are related to a few 
residual ambiguities of rm at 21 cm, and to the assumption 
that the A 2 -linear law is valid over the entire wavelength 
range. 

We co mpared the present RMs with those provided in 
Paper IVI . For 9 out of the 30 sources in that paper, the 
old Rotation Measure is rejected, while for 4 sources not 
detected at 13 cm, we can neither confirm nor disprove 
the old values. 

The source-frame Rotation Measure, (RMgr), is cal- 
culated multiplying the Galaxy-corrected RM □ by the 
factor (1 + z) 2 . 

6. Discussion 

6.1. The depolarising Faraday curtain: an empirical 
model 

To describe the Cotton Ef fect at the different observing 
wavelengths, in iPaper IVI . we outlined a model that re- 



3 We estimated the Mil ky Way contribution from the data set 
of Rotation Measures of iKlein et all (|2003l ), which contains 
about 200 sources carefully measured at 4 frequencies in the 
same sky area of our sources. From these data, the Galactic 
Rotation Measure is (—33 ± 7) radm -2 in the range 23 h < 
RA < h , (-80 ± 8) radm -2 in the range h < RA < 3 h , and 
(+9 ± 1) radm -2 in the range 7 h < RA < 15 h . 
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lates the source frame RM dispersion, <7RM, s f , to the source 
projected Linear Size LS, with a suggestion also of a de- 
pendence on z. 

More accurate values of <7Rjyr,sf> det ermined in the 
present paper, and further data from iPaper IV . are 
now available for 33 out of 38 of the B3-VLA CSS 
sources, which have Linear Sizes LS > 2.5 kpc, and 
known redshifts (either spectroscopic or photometric). 
Four of the missing sources (B3 0255+460, B3 1128+455, 
B3 1241+411, and B3 2349+410) are strongly depolarised 
at A > 6 cm or not polarised already at 3.6 cm, and their 
CRM.sf are likely to have a high value 1000 radm~ 2 ). 
For the last source (B3 1025+390B), ctrm, s i is not well 

A remaining sixteen 



constrained (see notes in Sect. IA. 1|) 
sources with LS > 2.5 kpc do not have a redshift and 
cannot be used in this analysis. We are however confident 
that their absence does not introduce bias in the results. 
Thus the sub-sample for which LS > 2.5 kpc, and known 
z, is adequate to investigate the dependence of <JRM,sf on 
both Linear Size (projected) and redshift. 

Sources smaller than 2.5 kpc are mostly depolarised 
at all frequencies. Therefore in the following analysis, we 
explore only the range LS > 2.5 kpc. 



We proceed as follows: 

a) The plot of CRM, s f versus LS (see Fig. |H1 middle panel) 
shows that: 

1- ORM.sf decreases with LS, with a large dispersion 
(about a factor 8 at » 90% level); 

2. there is a clear segregation in redshift, the high- 
redshift sources exhibiting, at the same LS, larger 
values of ciRM.sf ■ A hint of this dependence on red- 
shift is also seen in the distribution of ctrm.ods (see 
Fig. O top panel). 

To parametrize these findings, we assumed a 
model for the source frame RM dispersion that 
has a power-law dependence on LS and (1 + z), 
crm, mod <x LS~ a [l + z) b , and made a first estimate 
of a from the plot, neglecting at this stage the de- 
pendence on (1 + z). In this way, we obtained a 
first approximation model: CRM,modi <x LS . 

b) Using this law, we computed for each source a 
crm, modi [LS), corresponding to its Linear Size. These 
crm, modj take into account the dependence on LS, but 
not on redshift. The ratios <7RM, s f /cRM,modi show a 
large spread, part of which is due to the redshift de- 
pendence. A plot of these ratios versus redshift (see 
Fig. El bottom panel) allows us to estimate the pa- 
rameter b. We find b as 3.5. We note, however, that 
this power-law fit slightly underestimates CRM.sf at in- 
termediate redshifts, and overestimates CRM,sf at high 
redshifts. 

c) We used a second-order model, CRM,mod 2 ^ LS~ 2,5 (1 + 
z) 3 - 5 , computed for each source <7RM,mod 2 {LS, 1 + z), 
and analysed the ratios CRM,sf/cRM,mod 2 , as a function 
of both LS and (1 + z), and adjusted the two parame- 
ters a and b to derive another crm mod- 
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Fig. 9. Top: Plot of crm.ods vs LS. Middle: Plot of CRM,sf 
vs LS; the dotted line is the first order fit. Bottom: Plot 
of CRM,sf j corrected for the LS dependence, vs (1 + z); the 
dotted line is the first order fit. Different symbols are: (+) 
for z < 1.25, (A) for 1.25 < z < 1.8 and (A) for z > 1.8. 



d) We find a = 2.0 and b = 3.3. However, (see Fig. HOD the 
(1 + z) power-law fit still underestimates CRM, s f at in- 
termediate redshifts, and overestimates them at high 
redshifts: at intermediate redshifts 8 out of 10 data 
points are above the CRM,sf/cRM,mod = 1 line, while 



Rossetti et al.: B3-VLA CSSs: polarisation 



11 



at high redshifts 8 out of 10 are below this line. We 
evaluated the statistical significance of these system- 
atic effects by using both a contingency table and by 
analysing the medians and their variance in the two 
redshift bins. We find that the probability that the 
observed systematic effects are due to chance is < 5%. 
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Fig. 10. Plot of the ratios crRM,sf/c mo d vs (1 + z) for the 
power law model (1 + z) 33 . Symbols are as in Fig. [5] 



e) To remove this systematic effect, we split the sample 
into two redshift bins: (i) z < 1.8, and (ii) z > 1.8. 
We fitted the data for the lower redshift bin using the 
power-law relation (l + z) b , where b — 4.0; data for the 
higher redshift bin were fit using a constant, i.e. 6 = 
in the power-law relation. 

The adopted model is: 



f LS\~ 

ORM.mod = A ■ f — J (1 + z) 4 for Z < 1. 



\kpc 



(2.8) 4 for z > 1. 



(3) 



with A — 1.0 • 10 3 radm~ 2 kpc~ 2 . The uncertainties on a 
and b are ±0.25 and ±0.4, respectively. The parameters 
a and A may be slightly underestimated, because we did 
not account for the 4 small size depolarised sources, which 
probably have a high ORM.obs- 

In Figure 111) we compare the data and our adopted 
model from Eq. [3] we observe no systematic effects with 
respect to <7RM,sf / CRM,mod = 1, nor any significant depen- 
dence on LS or (1 + z). The residual dispersion of (7RM,sf 
about the model is a factor of approximately 2.5 at 90% 
level, and is likely due to intrinsic differences between one 
object and another. 



The model that we introduced in iPaper IVI was de- 
scribed by the following King law: 

-fold 



CRM, old 



[1 + {R/r c y\^ 



(4) 



RM.sf 
1 l~ 



J I I I III 

3 4 5 6 78910 



LS (kpc) 




1+Z 

Fig. 11. Top: Plot of the final model crRM,mod vs cjRM,sf • 
Middle: Plot of the ratios crRM,sf/cRM,mod y s LS. Bottom: 
Plot of the ratios ORM.sf/cRM.mod vs (1 + z). Symbols are 
as in Fig. \§\ In all three panels, the two dash-dotted lines 
are a factor of 2 above and below the model, and include 
> 80% of the sources. 



where R (kpc) is half the source linear size (R = LS/2). 
No dependence on (1 + z) was introduced, although 
the value of F id was assumed for the median redshift 



^median 



1. We derived the following parameter values: 
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-Fold = 1-5 x 10 5 radm -2 , 5 id = 2, and assumed the core 
radius r c — 0.5 kpc. 

For LS r c , Eq. 0] can be approximated by a power 
law similar to Eq. [3j 



CRM, old ~ -fold f c 4 {LS/2Y 



A \dLS 



(5) 



and for r c = 0.5 kpc, we find that .Add. — Fold- 

The two laws differ in their factor A and in their de- 
pendence on LS, the new law being much flatter. They 
intersect at LS « 3.8 kpc. 

The reason for the l arge differ ence between the old and 



new model is that, in iPaper IVI . the <7RM,oid(LS') model 



was derived using the Cotton Effect itself, namely from 
the (visually-estimated) A-dependent "critical LS"' below 
which the radio sources are almost totally depolarised. 
The steepness of the relation was strongly constrained by 
the "critical LS" at 21 cm. The new data at 13 cm have 
changed this situation. On the one hand, they have shown 
that, in a significant number of sources, a fraction of the 
polarised radiation remains constant at long wavelengths, 
and, on the other hand, they have increased the value of 
fRM,sf f° r a number of sources. The "critical sizes" are 
now the result of a combination of a dependence of CRM.sf 
on LS, and of partial coverage effects. Actually a num- 
ber of sources are still polarised at the longer wavelengths 
(A > 13 cm) as a result of a "partial coverage" (which 
becomes important for LS > 6 kpc), when the tTRM.sf, ex- 
pected for these sizes at the long wavelengths, according 
to Eq. [31 would be sufficiently large to totally depolarise 
all radiation for f c = 1. This is true, in particular, for 
21 cm. 
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6.2. The Cotton Effect Revisited 

We re-examined the "Cotton Effect" by constructing sev- 
eral model Cotton plots. These curves represent the de- 
polarisation, DP\ = m\/mo, that is expected at the ob- 
served wavelengths as a function of LS and z. 

We used different values of A, i.e. 500 < 
A (radm -2 kpc -2 ) < 2000, to take account of the dis- 
persion in 0RM,sf / crm, mod (see Eq. [3] and Fig.QT]). These 
choices provide values for cxRM,mod that are in the observed 
range of ORM,sf , for the majority (~ 80%) of sources. 

In Fig.rrjl we plotted only two family curves, i.e. those 
with z = 0.3, A = 500 radm -2 kpc -2 , and z = 1.8, A = 
2000 radm -2 kpc -2 , which represent a sort of minimum 
and maximum range for DP\ as a function of LS. 

The source depolarisation, DP\, lacking a measure of 
mo, is approximated by m\/m,3.6, with the assumption 
that TO3.6 is sufficiently close to the intrinsic degree of po- 
larisation. These data are plotted in Fig. [12] for two red- 
shift bins (z < 1.25, and z > 1.25). Sources without red- 
shift are also plotted. We label the polarisation as U DP\" 
to highlight that this is not the true depolarisation. This 
approximation, however, rises the following problems: 

a) When the source is undetected or barely detected at 
3.6 cm, it is undetected or barely detected at other 
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Fig. 12. From top to bottom: U DP\" at 6, 13, 21 cm (see 
text). In all panels (O) are for z < 1.25, (o) for z > 
1.25, (A) for Empty Fields. Sources showing an oscillating 
behavior are plotted as filled symbols. The upper curve is 
computed for z = 0.3 and .4=500 radm -2 kpc -2 in Eq.[31 
the lower curve for z — 1.8 and .4=2000 radm -2 kpc -2 . 



wavelengths, and "DP\ V would be the ratio of two 
small random numbers (noise/noise). This situation is 
particularly common for sources of small Linear Sizes. 
We believe that these sources are subject to strong de- 
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polarisation at all wavelengths, and their data points 
are plotted with the fictitious value U DP^ = —0.01, to 
emphasize that they are strongly depolarised. 

b) In the few sources of small Linear Size detected only at 
3.6 cm, the computed U DP\" could be overestimated 
because the measured TO3.6 could already be affected 
by depolarisation. 

c) We mentioned in Sect. 15.21 that some sources (per- 
haps as many as 9) display an oscillatory behaviour of 
m\, which we interpret as internal beating of at least 
two sub-components of different values of RM. In these 
sources a DP\" does not represent the depolarisation. 
However, for completeness, we plot these sources in 
Fig. [T31 using different symbols. 

The curves plotted in Fig. fT2l include the vast majority of 
data points at all wavelengths. Notable exceptions are: 

— a number of points significantly above the curves at 
small LS (typically 3-^5 kpc, depending on A). These 
sources, according to the model of Eq. [31 are expected 
to have a high <7RM,sf> and be totally depolarised, but 
this is not the case. We examined these objects one by 
one, and found that three with "DP 6 " significantly dif- 
ferent from zero have <7RM,sf smaller than average. At 
21 cm, the four sources with significantly high "DP21" 
have f c <^ 0.8, and therefore part of the polarised emis- 
sion escapes unpolarised. 

— Five points, typically at LS > 10 kpc, are well below 
the lower curve in Fig. \12\ as seen in the top panel 
( u DPq"). According to Eq.[3j these would be expected 
to have a small crRM.sf and hence be depolarised by a 
small amount, which is not the case. They may have 
an intrinsic value of ORM,sf that is much larger than 
average (A 3> 2000 radm -2 kpc -2 ). Only one source 
has a redshift and indeed it has the second highest 
CRM.sf m the sample. If this were the case also for 
the remaining four sources, these should disappear at 
other wavelengths. This effectively happens. The fact 
that these are mostly empty fields suggests that they 
might be at high redshift. 

— A few sources have "DPa"^ 1. These are sources with 
an oscillatory behaviour and m.3.6 depressed by the 
oscillations. 



6.3. The Faraday curtain: a physical model 

To provide a physical basis to the empirical model of 
the "Faraday curtain" (Eq. |3J), as in iPaper IVI . we adopt 
the model of a magneto-ionic medium (either smooth or 
clumpy) which is spherically symmetric about the radio 
source, described by a King-like profile of the relevant pa- 
rameters, and a randomly-oriented magnetic field. We as- 
sume, for simplicity, that the radio-source axis is orthog- 
onal to the line of sight (non-orthogonality effects will be 
considered in Sect. 16. 4j) . Any line of sight (los) to a point 
of the radio source will pass through several different ele- 
ments of the medium, of average size d much smaller than 
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Fig. 13. Sketch of the magneto-ionic ISM model around 
the sources. The line of sight (LOS) is the x axis. 



the source size, in each of which the (randomly-oriented) 
magnetic field B and the electron density n may differ 
from each others' (see a sketch of this model in Fig. [T3"|). 
Each medium element rotates the polarisation, according 
to the Faraday law, by an amount 5RM = knB\\d^\ As 
the field orientation changes at random from clement to 
element, the overall rotation along the line of sight (x co- 
ordinate) is on average zero, with a variance given by 



'RM.los 



d 2 n{xf B 2 \(x)dN, 



(6) 



where dAf = dx/d is the number of elements crossed by the 
line of sight in the interval dx. The integration is over all 
elements along that line of sight. For the smooth medium, 
this number is roughly the integration length divided by 
the element size. 

For both n(r) and B{r) we also assume a King-like 
distribution, i.e. 



n(r) 
B(r)-- 



n 



[1 + (r/rc) 2 ] 3 / 3 / 2 
Bp 

[1 + (r/r c ) 2 p/2 ' 



(7) 



where r c is the core radius of the medium distribution. In 
this model, orm.Ios is o btained by solving t he integral in 
Eq. O and, as shown bv lDolag et al. ( 2001 ). is given by: 



Csf,los(-R) = 



F 



Frf R 



(R/r c ) 

25 



2ii 



for R^> r c 



(8) 



4 In the following, we express the magnetic field B in fj,G, 
the electron density n in cm -3 , and every linear size (LS, R, 
d c , difL, r c ) in kpc. With these units, the value of the constant 
k in the Faraday law is 810 radm~ 2 /iG -1 cm 3 kpc -1 . 
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The dependence of F and S on the ISM parameters, varies 
depending on whether it is smooth or clumpy. These two 
components coexist, and we have to understand which 
plays the major role. 

6.3.1. Smooth medium 

We assume a "continuum" distribution of ISM elements 
("cells") of low density, with scale d — d c , central density 
no = n sm ,o and central intensity of the randomly-oriented 
magnetic field Bq = i? S m,o m Eq. [7] In this case, the pa- 
rameters of Eq. [5] are: 



F = K B smfi n smfi rl /2 d]J 2 radm 
. 609 + /i) - 1 



K is a factor dependent on 5 and, in our units, is in the 
range 300-400. 

In the case of equipartition between magnetic and ther- 
mal energy (isothermal case), or of a constant ratio be- 
tween the two, [i = (3/2. 

6.3.2. Clumpy medium 

This medium could be, for example, the Narrow Line 
Region (NLR), which consists of high-density (central 
value no = «nl,o) magnetized clouds, g?nl in size, with 
a randomly-oriented magnetic field of central intensity 
Bq = Bnl,0j and volume filling factor f(r). To ensure 
pressure equilibrium between the two components, we as- 
sume, in the first Eq. [7j that /3nl is equal to of the 
smooth medium in which the NL clouds are embedded. 
To evaluate the integral of Eq. [6l we have to model the 
distribution of the number per unit volume of the clouds, 
-Wnl ( r ) ■ We assume that the latter scales with r as a King- 
like law: 



AWr) 



[1 + (r/r c ) 2 p/ 2 



The cloud filling-factor f(r) is related to N^(r) by 
the relation: 



f(r) = 



N NL n di 



NL 



/o 



[1 + ( r / rc )2]3e/2 [ 1 + ( r / rc )2]3e/2' 



where /o is the filling factor at r = 0. The total number 
of clouds along a line of sight, at a distance R from the 
centre, Af(R), and the covering factor of these clouds, f^ h , 
are given by : 



Af(R) = H 



-/Vnl d NL r c 



H ■ (Jorc/dmf) 

[1 + (i?/r c ) 2 ] (3e - 1)/2 ~ [1 + (i?/r c ) 2 ] (3£ - 1)/2 



NL 



1 - e 



-AT 



where H is a factor dependent on e and has a typical value 
of a few units. 



F 

S 



The parameters of Eq. [H] are: 

v id 1/2 jl/2 j-1/2 j 

= K S NIj ,o n N L,o ^ c «NL Jo radl 
[609 + A*) + 3c] - 1 



(10) 



As in Eq. [9[ K is a factor dependent on 5, in the range 
300-400. 

Equations [T0l are similar to Eqs.El with cloud parameters 
in place of those of the smooth medium, apart for the 
factor fo and e. These terms account for the fact that 
along any integration length £ only a fraction / is filled 
by clouds, and that their number density decreases with 
increasing r. 



(9) 6.4. Physical vs empirical model 



If the polarised radiation were produced at the outer edges 
of the lobes only (hot spots), with the approximation R 3> 
r c , the second Eq. [8] could be related to Eq. 02 for R = 
LS/2, which is the length of each lobe on the assumption 
that the source is symmetric. We would derive 5 = 1 and 
F = A (2r c )- 2 (l + z) 4 . 

However, the polarisation does not always originate 
from hot spots only. We made some simulations that as- 
sumed a constant polarisation brightness, distributed over 
a fraction of the lobe length from the outer edges towards 
the centre, ranging in value between (hot spot only) and 
1 (entire lobe polarised). 

The simulations illustrate that we can describe the ef- 
fects of polarisation across a fraction of the source axis 
by introducing a factor h = R e g/(LS/2), where R e s is 
the "effective distance" between the source centre and 
the "centroid" of the polarised radiation, and LS/2 is 
the length of each lobe for a symmetric source. The 
range of h is, with good approximation, between h ~ 1.0 
(polarised radiation concentrated at the hot spots), and 
h w 0.5 (polarisation uniformly-distributed over the lobe). 
In other words, one should relate crRM,mod(£<S') (Eq. [3]) to 
cr RM ,lo S (#) (Eq. El) for R = R cB = h{LS/2). We should 
therefore have: 



F 



CRM, mod — 



[l + (h(LS/2)/r c ) 2 } s 
F{2r c f(h- LSy 25 



radm 

for LS > 2r c . (11) 



The location and extension of the polarised emission 
within the lobes change from source to source; this is likely 
to be one of the reasons for the residual dispersion of the 
data points with respect to Eq. [31 In Fig. [TTJ we expect 
that the lines (XRM.sf/cRM.mod = 1 correspond to the av- 
erage value (h) ~ 0.75, and sources with larger (smaller) 
h will lie above (below) the line. Comparing Eq. [HI with 
Eq.[3]we find: 



1 and F = A (2r c )- 2 {h) 2 {\ + z) 4 
w 140 r t T 2 (l + zf radm -2 . 



(12) 
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The condition LS 3> 2r c allows us to take r c = 0.5 kpc 
as a fiducial value for the core radius in what follows. 

We need to examine some further assumptions that we 
have made. In general, the sources are not in the plane of 
the sky and their shape may be asymmetric with respect 
to their cores, either intrinsically or because of an asym- 
metric ISM distribution. The orientation and asymmetries 
of the source structure may also cause the residual dis- 
persion in the data with respect to Eq. [3J 

We simulated different source orientations and arm- 
length asymmetries using Eq. [3] as an approximation to 
Eq. [U For sources that are not perpendicular to the line 
of sight, we found that the larger depolarisation of the far 
component is compensated by the lower depolarisation of 
the nearby component; this implies that the global <7RM,sf 
does not differ significantly (generally within 15%) from 
that expected for a source of similar LS, but in the plane 
of the sky. There are, however, some systematic effects 
that occur at orientation angles > 40° in sources of large 
sizes and/or at high redshifts. One consequence is that 
the parameter 8 could be underestimated and be > 1.1. 
A second consequence concerns the interpretation of the 
"partial coverage" and will be discussed in Sect. 16.81 

The arm asymmetry does not appear to cause major 
effects. 

6.5. Towards smaller sizes and higher frequencies 

It is interesting to use the depolarisation model of Sect. 16.31 
to predict what could be found at source sizes < 2.5 kpc 
and/or at wavelengths < 3.6 cm. 

The model shows that at short wavelengths the 
Faraday curtain may stop being foggy at small LS. The 
important parameter in this regime is the core radius, r c , 
for which we have adopted a tentative value of 0.5 kpc at 
large LS (Sect.[0]). 

Using Eqs. 1111 and Eq. [12] to compute (7R,M,mod, with 
r c = 1 kpc, we find that at 3.6 cm, for z < 1 and h = 0.5, 
DP never falls below 0.4 for LS > 1 kpc, and DP > 0.35, 
even for LS < 0.4 kpc. This contradicts the low level of 
polarisation found in our sources at LS < 1 kpc, which, 
instead, is compatible with the value r c < 0.5 that we 
have assumed. 

At 2 cm, from our model with r c ~ 0.5 kpc, we would 
expect DP > 0.2 at all sizes for z < 1. To our knowl- 
edge, integrated fractional polarisation data are rare in 
the literature at this frequency. At most, peak polarisa- 



tion flux-densities are reported (e.g. IStanghellini et al. 



200lh . We infer nevertheless, from these sparse data, that 
an r c <C 0.5 kpc may be required. 

6.6. Physical Properties of the Faraday curtain 

From the comparison of the empirical model (Eq. [3]) with 
the physical model (Eq. [8]), we can estimate the physical 
parameters for either the smooth- or the clumpy-medium 
component, assuming each to be independently responsi- 



ble for the RM dispersion, ORM,sf- We attempt to deter- 
mine which medium is causing the depolarisation. 

• The Smooth Component model. 

From the comparison of Eq. [9] and Eq. [12] we derive: 



K B smj0 n sm ,o (r c d c ) l/2 



14 °( 1 2 + ^ 4 radm- 

' C 

6(/? + M )-l i 



(13) 



where K « 380. From the second Eq.[l3l we derive /3+/i = 
5/6. If we assume that [i — (3/2 (see Sect. 16. 3|) . we derive 
(3 5/9. If instead B is constant (p, = 0), then (3 = 5/6. 
Therefore, the density scales with distance as: 



.(0 = 



[1 + (r/rc) 2 ] - 83 - 1 - 25 



The parameter (3 is in the range of that derived, using X- 
ray observations, for the hot component of the interstellar 
medium of nearby early- type galaxies. 

From the first Eq. Q2] we have: 



,0 B s 



cm 3 /iG 



0.37(1 + 2) 



kpc 



-1/2 



-5/2 



kpc. 

Introducing a constant ratio, g sm , between thermal 
(2n sm fi Icb T) and magnetic (i3 2 m0 /87r) energy densities, 
we derive: 

!^£ = i. 4 . io-5, 



cm 



B, 



sm,0 



TV" 



where T-j is the temperature of the smooth medium in 
units of 10 7 K. Combining these two expressions, we ob- 
tain: 

5=4 « i.s.io-^i-m 8 / 3 f^y^V 73 

cm 3 Vkpc / 



B s 



30 ft 



'1/3 



( d c r 5 c 



\kpc 



-1/6 



T 7 1/3 (l + z) 4 / 3 . 



We have no information about d c , apart from that it 
should be sufficiently small, compared to the sizes of the 
smallest sources, to produce the observed strong depolari- 
sation; therefore presumably it is on parsec scales. Hence, 
for r c < 0.5 kpc, d c ~ 1 pc, and Tj < 1, we have 
n sm fl(z — 0) w OAgll? cm~ 3 . For g sm w 1 (magnetic 
and thermal energy close to equipartition) , n snl) o is in the 
range of the values derived from X-ray observations for 
the hot component of the interstellar medium in the cen- 
tral regions of nearby early-type galaxies. 

In addition, the dependence of er s f,mod on (1 + z) 
(Sect. 16. ip indicates that the density of the medium in- 
creases strongly with redshift and then saturates. This oc- 
curs close to the "magic" redshift at which radio sources 
reach their maximum space density and perhaps "star for- 
mation" is close to its maximum. 

• The NL region model. 
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In IPaper we assumed that the steep decline of 
0RM,sf was due only to the decrease in the number of 
clouds per unit volume, as a function of r; we found that 
the implied covering factor of the clouds becomes negli- 
gible at sizes greater than a few kpc (the NL covering 
problem), making the model unrealistic. Using the new pa- 
rameters from the revised empirical model, and assuming 
a King-like distribution for cloud density and magn etic- 



Paper IVI . we 



field strength (Eq. [7]), not considered in 
re-examine the case. 

From the comparison of Eq. [T^] and Eq. 1101 we derive 
the following relations: 



KB NLfi ™NL,0 (fc^NL fa. 



1/a 140(1 + *) 



radm 



5 _ [6(/? + M )+3e]-l ^ 1(11) 



(K ps 380). 

To keep the number of parameters as low as possible, 
we assume that e = (3 and fi — (3/2 (as assumed for the 
smooth medium); hence, from the second Eq. 1 141 we ob- 
tain (3 = 5/12. If we had assumed B to be constant, we 
would have derived 8 = [(6/3 + 3e) - l)]/4 and (3 = 5/9. 

We observe that the density of the ambient diffuse 
medium, in which clouds are embedded, decreases in a 
way that is similar to the "smooth medium only" model, 
for a slightly lower value of (3 (0.43 < (3 < 0.56). For 
these values of {3, the number of clouds along the line 
of sight, Af(R), decreases slowly with R (oc R" 1 ^ 4 ), out- 
side the core radius, and f^ h is almost constant. If we 
assum e, instead , that ti.nl and -Bnl are both constant, 
as in IPaper IVI . then Af(R) drops quickly as a function 
of R, and f™ L becomes negligible, thus causing the NL 
covering problem. 

We now evaluate the strength of the magnetic field. 
From the first of Eq. [TJ] we derive: 



-Bnl.o 



0.4. 



cm 



-3 



kpc 3 



(/or c /d NL )- 1/2 (l + z) 4 



We first consider the situation for z = 0. We take 
^nl,o( z = 0) pa 5 x 10 3 cm~ 3 , w hich is in the range of 



value s quoted in the literature (e.g Peterson 19971 : iKoski I 
1978t h The quantity (/ r c /d^) is the average number of 
NL clouds along the line of sight, within the core radius. 
To have a covering factor > 0.9 for LS < 4 kpc, as re- 
quired by the data (Fig. [5]), (/o^c/^nl) has to be < 1. 
Taking r c < 0.5 kpc and fo ~ 10~ 4 (a valu e generally 
quoted in the literature, e.g. |Petersonl[i~997r ). we derive 
c^nl i> 0.05 pc, of th e order of what quoted from spectro- 
scopic observations ([Peterson II1997T ). We obtain: 



B Nh , (z = 0) « 6 iiG. 

The ratio between thermal energy and magnetic energy 
is: 

Snl « 10 4 . 



This model is not affected by the NL covering problem 
we experienced in IPaper IVI . It does, however, have some 
implications: 

i) The density inside the clouds, «nl, decreases with r, 

and at 5 kpc is already a factor 10 below the central 
value. 

ii) The magnetic field is quite low in these dense clouds, 
because the magnetic energy is a negligible fraction 
(« 10" 4 ) of the thermal energy. 

iii) The quantity ^nl-Bnl ^ s a strong function of z. If <7nl 
were independent of z, tt.nl and -Bnl would increase 
with redshift in proportion to (l + z) 8 / 3 , and (1 + z) 4 / 3 
respectively, out to z = 1.8 (where ORM.sf saturates). 
At this redshift, the cloud density tinl would be larger 
by more than a factor 10. If instead the magnetic-field 
strength were independent of z, the increase of ?1nl up 
to z — 1.8 would be a factor ps 60. 

Whether or not these implications of the NL model are 
acceptable is not clear to us. 

6. 7. Implications for young radio source evolution 

CSS sources are considered to be mostly young radio 
sources with typical ages far smaller than 10 6 years. They 
nevertheless form a large fraction of sources in radio cat- 
alogues. This old problem can be overcome by assum- 
ing that in early life-stages, the GPS/CSS phase, radio 
sources are luminous and then dim, due to adiabatic 
expansion moderately balanced by a continuou s energy- 



injection from the "nuclear engine" (see, e.g., |Sc heuer 



1974 [ [Baldwin Hl982tlFanti et al. lll99alReadhead et al. 
1996: Bcgelman 111 



In these simple physical models the evolution of the 
source luminosity depends on the density distribution of 
the ambient medium, which is assumed to be a power law 
[n(r) oc r~ n , with r\ — 3/3]. The range of values required 



by th e models is 1.5 < r] < 2.0 (see, e.g.. iFanti fc Fanti 
20031 and references therein). 

In our analysis of depolarisation parameters for a 
"smooth medium" , the density distributions derived are 
consistent with our above estimates. The density distri- 
bution for the NL model is flatter (1.2 < r\ < 1.6), but 
perhaps in agreement, within the uncertainties. 



6.8. Origin of "Partial Coverage" 

iBurn I ( 19661 ) discussed a variant of his own model (Eq.[l} 
that considered "partial coverage" . He assumed that the 
Faraday depolarisation is due to discrete cloud^l whose 
average number along the line of sight is N ' ■ If M ^> 1, 
the depolarisation is similar to that of the Gaussian model 
(Eq. [T]). If J\f <C 1, a number of lines of sight, how- 
ever, will not intersect any cloud. Therefore a fraction of 
the source's polarised radiation will emerge undepolarised 



5 Burn supposed that the clouds are in our Galaxy, but the 
model is easily usable for a location around the radio source 
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0039 - A 12 



A 2 (m 2 ) 

Fig. 14. B3 0039+412: example of data fit with Burn's 
Eq. [15] (dashed line) and with our Eq. [2] (dotted line) . 



through "holes" in the "curtain" keeping a constant level 
of fractional polarisation at long wavelengths. The rele- 
vant equation is: 



to = moe 



-7V(l-e 



(15) 



where F c is the RM of a single cloud. 

This equation, although formally different from the 
empirical one that we have used (Eq. [2]), is similar in 
shape. The differences cannot be discerned using the avail- 
able data, because of the limited wavelength coverage and 
the accuracy with which to a are measured (Fig. [T4|) . 

The revised Burn model would appear appropriate if 
depolarisation is due to clouds of the NLR. However, as 
discussed in Sect. 16. 6[ if the NL model that we have de- 
scribed is tenable, we expect that N(R) slowly decreases 
with R, and we expect covering factors > 0.9, while we 
also find values as low as 0.6-0.5 or less. 

Therefore we have considered other possible interpre- 
tations of the "partial coverage" , namely orientation ef- 
fects, intrinsic asymmetries in the radio-source structure 
and/or asymmetries in the distribution of the ambient 
medium. 

A promising alternative is related to effects of source 
orientation with respect to the line of sight. When the 
source axis is not perpendicular to the line of sight the two 
lobes suffer different depolarisations. In small-size sources 
of typical size < 4 kpc, the two lobes are rapidly depo- 
larised at A <C 13 cm, for any orientation and redshift. The 
overall depolarisation corresponds to the average 0RM,mod 
of the two lobes, which is close (within ~ 10%) to that for 
the case of orthogonality to the line of sight. For large 
source-sizes (> 8 kpc), deviations > 30° from the plane 
of the sky, and moderately high redshifts, cause the far 
component to be depolarised at A < 13 cm, and the near 
component little depolarised even at 21 cm. The effect in- 
creases with increasing LS, inclination angle to the plane 
of the sky, and z. The overall behaviour of depolarisation 
is a drop at short wavelengths, followed by a much slower 




Fig. 15. Simulation of a symmetric radio source at 45° 
from the plane of the sky with LS = 8 kpc and z = 0.6. 
The two Gaussians with too = 0.5 (dashed and dotted 
lines) represent the depolarisation of the far and near lobes 
(respectively strongly and weakly depolarised) as a func- 
tion of A 2 . The dash-dotted line is the sum of the two 
contributions. The "experimental" data points are values 
read from the model curve at the observed wavelengths 
of the present experiment. The continuous line is the fit 
to the four data points using Eq. [3] with f c = 0.78, i.e. 
assuming a partial source coverage. 



decrease (see e.g. dash-dotted line in Fig. [15]) that, when 
allowance is made for the errors, is indistinguishable from 
a flattening. We made simulations for h = 0.75 in Eq.fTTl 
and found that, in a number of situations, the depolari- 
sation, with only the four wavelengths available to us, is 
well fitted by Eq.H with 0.5 < f c < 0.8. 

We analyzed also the possible effects of intrinsic asym- 
metries in the source arm length. The two lobes suffer 
different depolarisations because the shorter arm, located 
closer to the centre, is inside the inner denser region of 
the medium (see Eq. [TT]). We made simulations base d 
on the arm-ratio distribution of Rossetti et al. I (|2006l) . 
for the 21 B3-VLA CSS sources for which a core is de- 
tected at 15 GHz. We found that, under the assumption 



of a spherically-symmetric Faraday medium, no significant 
effects are expected. The results might be different if the 
source asy mmetries were caused by density asymmetries 
in the ISM (jjevakumar et al. 1120051 ) , the denser ISM pre- 
venting one of the two lobes to grow like the other. In this 
case, the cTRM.sf of the two lobes would differ not only 
because the shorter lobe is closer to the source centre (as 
before) but because, in addition, the ISM is denser on the 
short lobe side. We did not model this situation because it 
requires too many parameters that cannot be constrained 
easily. We used, instead, an experiment al approach. For 
13 of 21 sources observed at 15 GHz ( Rossetti et al. I 
20061 ) that belong to the WSRT sub-sample, we plotted f c 
versus arm ratio. No significant relation has been found, 
suggesting that, if radio source asymmetries are due to 
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ISM asymmetries, the latter are not very important in 
simulating a partial coverage effect. 

In conclusion, we suggest that covering factors smaller 
than 0.8, are possibly justified by orientation effects. 

6.9. Rotation Measures: A large-scale ordered 
magnetic field? 

In the top panel of Fig. 1161 we plot <7RM,obs versus the ob- 
served, Galaxy-corrected, Rotation Measures, RM ^ S for 
the 30 sources (21 with z) for which both data are avail- 
ably. We discuss possible selection effects in the plot. 
Sources for which URM.obs > 70 radm" 2 and f c > 0.9 
should be strongly depolarised at 13 cm (Eq.[JJ), therefore 
their RMs are not measurable because the detectability 
at this wavelength is the condition we adopted to com- 
pute RM (Sect. [5~3j) . The 10 sources in the plot that have 
ORM.obs > 70 radm -2 , have f c < 0.9. 

Because of band-width depolarisation, most sources 
with RM ohs > 750 rad m -2 would be depolarised at 13 cm 
and their RM would not be measurable. The vertical and 
horizontal lines show these limits. In the bottom panel of 
Fig. 1161 we plot similar data, with the source frame pa- 
rameters. 

Both figures show that about half of the sources have 
RM > ctrm- This is unexpected for a model in which the 
magnetic field of the Faraday screen is random on small 
scales. Our Monte Carlo simulation (Sect. I5T2")) shows that 
in these models the E vector is expected to have, statis- 
tically, a small global rotation, with \RM\ < orm/V^c 
(where N c is the ratio between source area and cell area) 
as long as ct rm A 2 < 2 rad (DP p> 5%). If N c > 100, the 
data points in Fig. [TBI would be expected to cluster around 
or above the uppermost line, and, given the values we find 
for ctrm (see Table [A~Tj) . no large rotations of the polar- 
isation angle with A 2 , the residual of the RM dispersion, 
are expected. In other words, with the possible exception 
of a minority of sources (< 15%), the RMs are not the 
residuals of the random rotations that depolarise the ra- 
diation. They are likely generated on a larger scale by an 
ordered magnetic field component. 

In Paper ivl . we found a marginal correlation between 
RM S [ and z. In the top panel of Fig. [T7l we plot RM s f as 
a function of (1 + z). For RM s f < 100 radm" 2 , the dis- 
tribution in z is uniform, while for RM s f > 100 radm" 2 , 
8 sources out of 10 have z > 1. There is a hint of an 
upper bound that increases with (1 + z). The probabil- 
ity that the observed distribution is a random selection 
from a uniform redshift distribution is ~ 7%. We fur- 
ther note that sources for which RM s f > (TRM,sf (circles 
in Fig. 1171 top panel) are distributed in a broad band 
such that RM s f oc (1 + z) 3 . The probability that this is 
caused by a random selection from an uniform distribu- 
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6 We recall that 15 other sources (10 with z) with derived 
crm are unpolarised at 13 cm so that their RM is unknown 
(see Sect. 15.3]) . In addition, 2 sources (1 with z) have 7?Mbut 
not orm 



RM S , 

Fig. 16. Plot of ctrm vs RM in the observed (top panel) 
and in the source (bottom panel) frame. Open symbols 
in the top panel indicate Empty Fields. The tilted lines 
represent (CTR M ,obs/-RM ob s) and (CT RM , s f/i?A/ s f) = 10, 3, 
1, 0.3, 0.1. 

tion is < 1.5%. We have to examine if there are any biases 
which can affect the source distribution in the plot RM 
vs (1 + 2;). There are 10 sources for which we do not have 
a measured RM because they are depolarised at 13 cm. 
If these sources had an RM above the apparent upper 
bound in Fig. |T7l (top panel), their RM would be high 
(> 300 radm -2 ) in the observer's frame; we would expect 
therefore to see a rotation of the polarisation angle be- 
tween 3.6 and 6 cm. At first sight there appear to be no 
candidates for these "missing RM' among the unpolarised 
sources at 13 cm. Although this statement has to be taken 
with some caution, we conclude that there is a real corre- 
lation between RM S { and (1 + z), at least for sources for 
which RM s f > CTRM,sf- 

In the bottom panel of Fig.[T7l we plot RM S { versus LS. 
The distribution of data points is peculiar. There seems to 
be a "sequence" of sources, for RM s f > 100 rad m" 2 , for 
which RM s f oc LS 1 " 2 2 , which is a dependence on LS that 
is similar to that of <7RM, s f ■ This "sequence" (9 sources out 
of 21) is composed mainly of objects for which RM s f > 
CRM,sf (8 out of 9) and z > 1 (7 out of 9). Of the remaining 
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Fig. 17. Plot of RM Bf vs (1+z) (top panel) and of RM sf 
vs £5 (bottom panel). In both panels (A) and (o) represent 
objects with RM s f < and respectively > cr s f ; filled symbols 
in the bottom panel represent objects with z > 1. 



objects, RM s f does not appear to show correlation with 
LS. 

We conclude that for 50% of the CSS (mainly at high 
redshifts, (z) w 1.5), a large-scale ordered magnetic field 
is present and produces the higher values of RMs. Its scale 
must be approximately 10-20 kpc to justify the correlation 
observed in the bottom panel of Fig. 1171 

Applying a model similar to that for (jRM.sf fSect. l(jT3l 
to the data of the "sequence" RM S { versus LS, we find 
that the required ambient density and magnetic-field de- 
pendence on r, are similar to that for ORM.sf (|3 + /i« 1), 



and B 



40 r: 



10 fj,G, at the median z « 1.9 



of the sources in the "sequence" . We observe that not all 
high-redshift objects lie on the RM a { oc LS~ 2 2 relation. A 
possible explanation is that in a large-scale magnetic-field, 
the orientation effects are important for RM. High-redshift 
sources that do not belong to the sequence, should have a 
line of sight that is at significant angles to the magnetic 
field. 

If Bnono oc (1 + z) 4 , as for ORM,sf, at z < 1 the RM 
would drop to < 150 radm -2 . Most of the low-redshift 



sources would fail to show the correlation of the high- 
redshift ones. 



7. Summary and Conclusions 

We have observed 65 radio sources from the B3-VLA CSS 



sample (jPaper J ) using the WSRT at 13 cm, to study the 
source polarisation properties. At the WSRT resolution, 
the radio sources are all unresolved, and we can therefore 
only discuss their global properties. The new 13-cm data, 
combined with earlier low-resolution VLA data at 3.6, 6, 
and 21 cm, have improved the determination of Rotation 
Measures (RM) and the Rotation Measure Dispersions 
(crm)- This has allowed their properties to be defined 
more carefully as a function of both redshift and Linear 
Size, and the characteristics of the surrounding ISM to be 
modelled. 

The main results of the paper are the following: 

1) Radio sources < 5 kpc are mostly unpolarised at 
13 cm. This result is similar to results fou nd at 20 cm 
( Cotton et al. Il2003h and at 3.6 and 6 cm ([Paper IVI ). 
and the "critical size" that we find is intermediate be- 
tween those found earlier. 

2) The 13 cm polari sation ang les have led to a revision of 
the earlier RMs ([Paper IVI ) for « 30% of the sources. 

3) The integrated fractional polarisation m\, as a func- 
tion of wavelength shows in general a decrease be- 
tween 3.6 and 13 cm, followed by a flattening between 
13 and 20 cm. This is the major observational result 
of the p resent pap er. At variance with the conclu- 



sions of IPaper IVl the iBurn I l| 19661 ) and iTribble 



( 19911 ) models are not adequate to describe t his be 



Burn 



haviou r. However, an empirical variant of the 
(1966) model, which introduces a "partial coverage" 
by the depolarising curtain, appears to reproduce the 
data well. The adopted formula allows the Rotation 
Measure Dispersion, crm, an d the fraction, / c , of 
source covered by the depolarising curtain, to be de- 
termined. 

4) For a minority of sources, m\ shows an irregular, pos- 
sibly oscillatory, behaviour with A. We propose that 
these sources contain sub-components with different 
Rotation Measures that produce beats with A in the 
integrated fractional polarisation. 

5) orm shows a clear dependence on redshift (up toz« 
1.8), and on projected Linear Size. 

6) We have analysed these de pendences using a model 
similar to that presented in IPaper IVI . in which a de- 
polarising Faraday curtain, with a King-like distribu- 
tion of the magneto-ionic medium, is produced cither 
by a smooth medium with an irregular magnetic field 
on small scales, or by a clumpy magnetised medium 
(NL region). The parameters der ived in eit her case 



differ from those we obtained in IPaper IVl . because 
of the use of the new 13 cm data and of the new model 
adopted for the depolarisation behaviour with A. In 
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both models the core radius of the curtain has to be 
< 0.5 kpc. 

In the smooth medium model, the required central 
density, at z = 0, is in the range of that found using 
X-ray observations in the centre of early-type galax- 
ies. Beyond the core radius, it decreases according to 
r~ v , where 1.5 < r/ < 2.0. The central magnetic field, 
at z = 0, is ~ 150 fiG, and the magnetic field energy 
density is close to the thermal energy. 
In the clumpy medium model, identified with the 
NLR, we assumed pressure equilibrum between the 
clouds and the diffuse medium in which they are em- 
bedded. The required magnetic field, inside the clouds 
at z — 0, is 6 fiG, and the magnetic energy is a negli- 
gible (w 10~ 4 ) fraction of the thermal energy. Beyond 
the core radius, the density of the smooth medium, 
which confines the clouds, should decline as r _7? , with 
1.2 < r\ < 1.6. The product -Bnl,o x «-nl,o is a strong 
function (oc (1 + z) 4 ) of redshift up to z » 1.8. 
The NL model requires more parameters than the 
smooth medium model, and these parameters need tun- 
ing to avoid the "NL covering problem" discussed in 
Sect. 16.61 For these reasons, we favour the smooth 
medium model. 

7) The ISM density profiles required by the depolarisation 
models are consistent with the expectations of the evo- 
lutionary models of young radio sources. 

8) The partial coverage, that we introduced to describe 
the long wavelength behaviour of m\ vs A 2 , is proba- 
bly due to orientations effects of the source with re- 
spect to the line of sight. For sources that are almost 
orthogonal to the line of sight, the two lobes are depo- 
larised similarly by a spherically-symmetric medium, 
and m\ is fitted by the Burn model. However, when 
the source is at large angles with respect to the line of 
sight, the far lobe is more depolarised causing the rapid 
decrease of m\ at short wavelengths, while the nearby 
one is less depolarised giving rise to a lower decrement 
of m\ with wavelength. The total fractional polarisa- 
tion m\ then appears to be almost constant at long 
wavelengths. 

9) The total source-frame Rotation Measures are gener- 
ally too large to be the residual of random rotations 
across the source due to a totally irregular magnetic 
field, and require a large-scale ordered magnetic-field 
component. The Rotation Measures show hints of a 
correlation with both redshift and Linear Size in a way 
similar to the Rotation Measure dispersions. If these 
correlations are real, the ordered magnetic-field com- 
ponent would have parameters similar to those of the 
random field component. 
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Name 


Sl3 




mi3 




X13 


RM obs 


RM st 


°"RM,obs 




m 


fa 


z 


LS 


Notes 


(i) 


(2) 




(3) 




(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


0034+444 


0.42 


1 


40+0. 


11 


-33 


-85 


-72 


41 


589 


3.1 


0.9 


2.79 


11.9 


PI, N 


0039+373 


0.60 





.10+0. 


10 








243 


982 


1.0 


0.85 


1.01 


0.5 


P5 


0039+398 


0.48 





.34+0. 


11 


42 
















10 


P6 


0039+412 


0.25 


2 


.81+0. 


12 


-2 


-100 




22 




5.2 


0.64 




8. 


5 


PI, N, m 


0041+425 


0.30 





28+0. 


13 


(-24) 






>368 




3.0 


0.93 




4 


.1 


PA 


0049+379 


0.46 





.27+0. 


11 


(17) 






254 


1852 


9.0 


0.98 


1.7 


6. 


1 


P3 


0110+401 


0.42 


3 


.54+0. 


11 


54 


-65 


-92 


~ 8 


~ 49 


5.4 


[1] 


1.48 


10 


.7 


PI, N, m 


0120+405 


0.38 


1. 


.41+0. 


11 


-102 


-340 


-880 


140 


474 


5.5 


0.75 


0.84 


10. 


1 


PI, N 


0123+402 


0.16 


1. 


.32+0. 


15 


-48 


-136 a 




<50 




4.5 


>0.95 




5. 


2 


P2 


0128+394 


0.20 





.55+0. 


14 


-87 


-150 


-473 


~ 10 


~ 68 


0.8 


[1] 


1.6 


11 


.5 


P5, N, m 


0137+401 


0.19 


1. 


.53+0. 


14 


■30 


-88 a 


-55 


<21 


144 


1.7 


> 0.4 


1.62 


17 





PI, N 


0144+432 


0.23 


1. 


.05+0. 


12 


-40 


-54 


+ 132 


105 


536 


3.5 


0.66 


1.26 


17. 


2 


PI, N 


0147+400 


0.51 





.01+0. 


10 





















4 


P6 


0213+412 


0.39 





63+0. 


11 


10 


-67 a 


+30 


63 


146 


3.1 


0.89 


0.52 


7. 


3 


PI 


0222+422 


0.16 





39+0. 


14 


(60) 






195 


3949 


4.0 


0.9 


3.5 


11. 


9 


PS 


0228+409A 


0.24 


1. 


60+0. 


12 


38 


-159 




246 




4.0 


0.65 




16. 


2 


PI 


0254+406 


0.33 


1. 


.53+0. 


12 


3 


-27 a 


+262 


<20 


<99 


1.7 


>0.25 


1.22 


18 


.9 


PI 


0255+460 


0.45 





.10+0. 


11 
















1.21 


3 


2 


NP 


0701+392 


0.35 





.15+0. 


12 








102 


512 


1.0 


0.8 


1.24 


7. 


.8 


PS 


0722+393A 


0.64 





.10+0. 


10 


















1. 





NP 


0744+464 


0.32 





.50+0. 


11 


-3 


-277 a 


-4417 


229 


3537 


5.5 


0.9 


2.93 


5. 


1 


PI, N 


0748+413B 


0.13 


1. 


.40+0. 


13 


13 


+29 




201 




18 


0.94 




1. 


.8 


PI 


0754+396 


0.33 


8. 


.32+0. 


13 


-97 


+5 a 


-39 


63 


613 


11.6 


0.32 


2.12 


8. 


9 


PI 


0800+472 


0.60 





33+0. 


10 


-54 


-19 


-64 


148 


337 


2.4 


0.83 


0.51 


3. 


.7 


PI 


0805+406 


0.37 


5. 


.75+0. 


12 


2 


+22 a 




30 




6.5 


0.3 




10. 


7 


PI 


0809+404 


0.75 





.15+0. 


10 








210 


519 


4.5 


0.97 


0.55 


4. 


5 


PS 


0810+460B 


0.66 





.11+0. 


10 








>270 


>478 


2.0 


0.94 


0.33 


3. 





PA 


0814+441 


0.17 


1. 


93+0. 


12 


-64 


+352 














17 





P5, N 


0822+394 


0.77 





.02+0. 


10 
















1.18 





2 


NP 


0840+424A 


1.03 





06+0. 


10 





















.3 


P6 


0902+416 


0.35 


1. 


90+0. 


11 


-65 
















1. 


4 


P5, N 


0930+389 


0.17 


2 


.79+0. 


12 


-45 


+46 


+428 


<24 


<277 


2.8 


> 0.8 


2.4 


14 


.4 


PI, N 


0951+422 


0.29 





.44+0. 


11 


-3 


+7 


-15 


206 


1592 


4.5 


0.88 


1.78 


8. 


3 


PI 


0955+390 


0.31 


6 


.27+0. 


13 


15 


+ 15 a 




<13 




6.6 


> 0.6 




10 


2 


PI 


1007+422 


0.29 





.13+0. 


13 








400 




6.5 


>0.95 




0.55 


P5 


1025+390B 


0.49 


2. 


98+0. 


11 


-101 


20 


-54 










0.36 


9 


.7 


P2, N, m 


1027+392 


0.28 





26+0. 


11 


(91) 



















9 


P5, N 


1039+424 


0.17 


6 


90+0. 


15 


-57 


+ 12 a 




<13 




7.3 


> 0.8 







.5 


PI 


1044+404A 


0.25 





35+0. 


11 


14 


+376 


(9546) 


155 


(+4032) 


8.5 


0.95 


4.1 


3. 


2 


P2, N 


1055+404A 


0.25 


5. 


69+0. 


12 





+8.9 a 




<13 




5.9 


> 0.6 




12 


.1 


PI 


1128+455 


1.35 





00+0. 


10 
















0.40 


2. 


9 


NP 


1136+420 


0.32 





.20+0. 


11 








>324 


>1085 


2.4 


0.9 


0.83 


4 


3 


PA 


1201+394 


0.32 


2 


29+0. 


11 


26 


+38 


+61 


<24 


<51 


3.0 


> 0.90 


0.45 


7. 


2 


PI, N, m 


1204+401 


0.15 





.31+0. 


15 


(-1) 






68 


641 


5.2 


>0.95 


2.07 


6. 


.7 


PS 


1216+402 


0.23 


2 


.03+0. 


12 


64 


-190 


-616 


<10 


< 30 


4.0 


[1] 


0.76 


15 


.5 


PI, N, m 


1220+408B 


0.28 


1. 


.47+0. 


11 


-100 


+28 a 




61 




6.9 


0.95 




17. 





PI 


1225+442 


0.24 





09+0. 


12 
















0.22 





.9 


P6 


1233+418 


0.49 


1. 


.15+0. 


11 


-52 


-45 


-84 


50 


78 


4.3 


>0.9 


0.25 


5. 


3 


P2, N, m 


1242+410 


1.04 





.04+0. 


10 
















0.81 





3 


NP 


1340+439 


0.32 





23+0. 


11 


(-21) 



















.3 


P6, N 


1343+386 


0.66 





38+0. 


11 


-15 


-379 


-3129 


115 


928 


3.8 


0.86 


1.84 


0. 


5 


PI, N, m 


1350+432 


0.09 


3 


.18+0. 


16 


-86 


+230 


+2193 


150 


1490 


10.0 


0.55 


2.15 


0. 


5 


PI, N, m 


1432+428B 


0.66 





.15+0. 


10 





















2 


P6 


1441+409 


0.67 





10+0. 


10 








>280 




1.4 


>0.9 







4 


PA 


1449+421 


0.43 





.12+0. 


10 





















.3 


NP 


1458+433 


0.29 


2 


.70+0. 


11 


12 


—5 a 


-52 


<15 


<55 


2.7 


> 0.2 


0.93 


6 


9 


PI, N 


2301+443 


0.69 





.03+0. 


10 
















1.7 


2 


.1 


P6 


2302+402 


0.83 





.10+0. 


10 


















2 


6 


P6, N 


2304+377 


1.07 





.19+0. 


10 








182 


357 


3.2 


0.96 


0.4 


0. 


3 


P5 


2311+469 


1.38 


4 


22+0. 


11 


-70 


-8 


+77 


41 


126 


5.5 


0.36 


0.75 


9. 


6 


PI, N 


2322+403 


0.22 





03+0. 


12 








207 




6.5 


0.98 




14 





PS 


2330+402 


0.60 





06+0. 


10 





















.3 


PC, 


2348+450 


0.50 





.04+0. 


11 








>290 


1137 


2.5 


0.95 


0.98 


1. 


3 


PA 


2349+410 


0.28 





.24+0. 


11 


(76) 














2.05 


4 


.7 


P5, N 


2358+406 


0.99 





02+0. 


11 


















0. 


3 


NP 



- Columns 1, 2 and 3: Source Name, 13 cm flux density (Jy), fractional polarisation and error (%); 

- Column 4: 13 cm Electric vector p. a. (degree) for sources with mi3 > 3a m ; values in parenthesis are for sources with 
2<r m < mi3 < 3<r m . Errors, computed as in Sect. 13. ll are < 9° for 772,13 > 3er m 

- Columns 5 and 6: Observed and source frame Rotation Measure (radm -2 ). Formal errors of RM h B are < 10 radm -2 . An "a" 
near RM b B means a good A 2 -linear fit (see Sect. 15. 3[ ). 

- Columns 7 and 8: Observed and source frame orm (radm -2 ); formal errors are typically < 20% 

- Columns 9 and 10: Intrinsic Fractional polarisation (%) and Covering factor 

- Columns 11, 12 and 13: Redshift (photometric one decimal digit only), Largest (projected) Linear Size (kpc), Notes 

PI: the source is detected at > 3ap at all four wavelengths; P2: the source is detected at > 3op at 3.6, 6, and 13 cm; P3: the 
source is detected at > 3op at 3.6 and 6 cm; P4: the source is detected at > 3<rp at 3.6 cm only; P5: the source is detected 
at > 3<rp at two or three non contiguous wavelengths; P6: the source is detected at > 3crp at 21 cm only; NP: the source is 



22 



Rossetti et al.: B3-VLA CSSs: polarisation 



Appendix A: Polarisation Data 

A.l. Notes to individual sources 



0034+444: See Fig.[5]and Sect. [El The x 3 . 6 is discrepant 
from the x(^ 2 ) behaviour at the other frequencies, causing 
a poor overall fit. We have not been able to improve the 
fit using the "two polarised component model" of [B] 
0039+412: The RM obs obtained by a AMinear fit, using 
also the individual IFs of the f 3 cm band, does not prop- 
erly account for X2i- We have used the "two polarised 
component model" , constrained by the polarisation struc- 
ture seen by the VLA at 3.6 and 6 cm, which introduces 
a small modulation in the polarisation angles over the A 2 
behaviour which better fits all the data. The RM D b s in 
Table [A~T1 is from this model. 

0110+401: The source shows a deep minimum in frac- 
tional polaris ation at ~ 11 cm (d ata from the Effelsberg 



telescope, see lKlein et al. 



(|2003l )V A AMinear fit, using 
also the individual IFs of the 13 cm band, does not ac- 
count for \2i- We have used the "two polarised compo- 
nent model" which fits very well both behaviours of m\ 
and xM vs ^ 2 ( see Fig. [6|). The RM b s and <7RM,obs in 
Table I A. II are from this model. 

0120+405: An alternative value of RM ohs is « —154 rad 
m~ 2 . The adopted RM obs justifies the lack of polarisa- 
tion at 21 cm as a bandwidth depolarisation effect. Both 
the reported RM obs and the alternative one are in bad 
agreement with RM13. 

0128+394: This source shows a possible oscillatory be- 
haviour of m vs A 2 . We have used the "two polarised com- 
ponent model" , constrained by the polarisation structure 
seen by the VLA at 3.6 and 6 cm, which fits both m and 
X vs A 2 . The i?M b s and <7R,M,obs in Table |A~T1 are from 
this model. 

0137+401: The angle X3.6 is in poor agreement with a 
A 2 -linear fit. 

0144+432: The fitted i?A/ bs has a poor chi-square, 
mainly due to X2i- No acceptable two component model 
has been obtained to improve the fit. 
0744+464: According to our RM obs the fractional polar- 
isation at 21 cm is depressed by a large factor because of 
bandwidth depolarisation. 

0814+441: The source is undetected at 6 cm and well 
detected at the other three wavelengths, suggesting a pos- 
sible oscillatory behaviour of m vs A 2 . The two polarised 
component model quite fits the data except X21 well. The 
fitted i?M b s , which is very large, would imply virtually 
no polarisation at 21 cm because of bandwidth depolari- 
sation, contrary to what is observed. However, the model 
shows strong modulations of i?Af b s , and at 21 cm the 
"local -RAiobs" is lower than the average, implying only a 
reduction of mix of a factor ss 2, in agreement with the 
data. 

0902+416: This is another source with a possible m\ os- 
cillatory behaviour, being well detected at 3.6 and 13 cm 
but not at the other two wavelengths. We have modelled 
it with the two polarised component model; however, be- 



cause of the small amount of data and lack of informa- 
tion on polarisation sub-structure, we consider the results 
too uncertain. No polarisation parameters are reported in 
Table Ell 

0930+389: The RM ohs quoted in Table ED is strongly 
constrained by RM13. The angle X3.6 is in bad agreement 
with the overall fit and causes a poor chi-square. 
1025+390B: No data is available at 21 cm. The fractional 
polarisation TO3.g is a factor w 2 lower than at 6 and 13 cm. 
It may be another case of oscillatory behaviour of m\. We 
have modelled the source with a two polarised component 
model, but, due to the limited amount of data, only RM obs 
is given in Table IA.ll 

1027+393: The source is strongly polarised at 6 and 
21 cm but not at the other wavelengths. Therefore it is 
another case of oscillatory behaviour of m\ . Although the 
two polarised component model fits both m\ and xWi 
because of the limited amount of data the results are not 
properly constrained. 

1044+404A: The very large RAl D b s is likely the cause of 
the lack of polarisation at 21 cm because of bandwidth 
depolarisation. 

1201+394: An improved fit of x(A) vs A 2 is obtained by 
the "two polarised component model" which introduces a 
small modulation of the xW behaviour. Also the fit of the 
amplitudes is improved. The data in Table I A. II are from 
this model. 

1216+402: A possible oscillatory behaviour of m\ is well 
described by the "two polarised component model" which 
well fits tji\ and x(^)i except for X2i- RM13 is somewhat 
lower than the adopted i?M b s , but the model justifies 
it as a result of modulation of xW because of the beat 
between the two components. Data in Table [Aj] are from 
this model. 

1233+418: The "two polarised component model", based 
also on the polarisation substructure seen by the VLA at 
3.6 and 21 cm, gives an improved fit to both m\ and xW- 
The parameters in Table [A~l1 are from this model. 
1340+439: This source is well detected at 21 cm and 
possibly detected at > 2ap at 3.6 cm. It may be another 
example of the oscillatory behaviour of m\ , which is easily 
fitted by the "two polarised component model" . However, 
because of to the poorness of the data we do not report 
any polarisation parameter. 

1343+386: The angle Ya.fi is very discre pant with respect 
to the overall fit. lOrienti et al. I (|2004l) showed that the 
polarisation is from the bright southern hot spot, and has 
substructures, with i?M Q b s differences ~ 100 rad m~ 2 . We 
have applied the "two polarised component model" and 
obtained a much better fit for m\ and xityi except for 
X3.6 which we consider an outlier. The data in Table |A~T1 
are from this model. 

1350+432: At a first sight there is a large discrepancy 
between RM13 and the global RM b s - However, the "two 
polarised component model" properly fits the data, and 
justifies the low value of RM13 as a result of modulation 
because of the beating of the two components. Data in 
Table IA.1I are from this model. 
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1458+433: m 3 6 = 1.8% is low compared to m 6 = 2.6%, 
m 13 = 2.7% and m<z\ = 2.4%. We have no obvious expla- 
nation for this behaviour and have htted ORM.obs with the 
exclusion of m 3 .g. 

2302+402: This source is detected at 21 cm only with a 
S/N = 17. The fractional polarisation is much larger than 
the upper limits at the other wavelengths. 
2311+469: the angle \3.6 is quite discrepant from the 
global X (A) fit. 

2349+410: Undetected in polarisation at 6 and 13 cm. 
Another possible case of oscillatory behaviour of m\, 
which can be fitted by the "two polarised component 
model" . However, because of the poorness of the data we 
do not report any polarisation parameter. 

Appendix B: The two polarised component model 

The polarisation status of a source can be represented by 
a vector of amplitude Sp and position angle 2\, defined 
in Sect. 13.21 When two sub-structures are present, differ- 
ent RMs in the two components may cause a differential 
rotation of the two polarisation vectors as a function of 
A 2 . The result is that the total source polarisation and E 
vector position angle, under vectorial addition of the two 
component polarisation vectors, may present signs of in- 
terference. This produces, for instance, minima in >Sp(A 2 ) 
when the two polarisation vectors happen to be opposite. 

We have applied this two polarised component model 
to a number of sources having one of the following char- 
acteristics: 

i) discrepancies between the m\ and the fitted model 

(Eq. [2]) , mostly possible oscillations of m\ vs A 2 , which 
could be indications of a beat between polarised sub- 
components with different RM; 

ii) poor fit of the polarisation angles with the A 2 -linear 
law, again a possible indication of modulation which 
are due to polarised subcomponents with different RM. 
These sources should also exhibit beats, but there is 
only a marginal evidence for them. However, the os- 
cillations might have been lost because of the limited 
A 2 sampling or might have strongly been damped by 
depolarisation. 

The parameters involved in the model are: the intrinsic 
fractional polarisation mo, intrinsic position angle \o and 
RM of each component, and a common <trm and f c which 
accounts for the depolarisation. The available data are the 
m\ and x(A), (typically 2x4 data), and RM 13 fSect. l3~2~4l 
when available, so that we mostly have a comparable num- 
ber of parameters and data points. It makes no sense to 
use more than two components in the model because the 
number of parameters would exceed the number of data. 

In order to constrain the model as much as possible 
we used as starting parameters, when available, the data 
of the two more p olarised co mponents from the high reso- 
lution VLA data (jPaper IVT l: i) m 3 . 6 , ii) i?M| 6 , Hi) X3.6 



Furthermore it is also possible that the beat is due to sub- 
components unresolved even at our best VLA resolution, 
so that the input data would not be the relevant ones. 

Sources of group i) are mostly fitted by the model, 
ranging from good and well constrained cases (see e.g. 
B3 0110+401) to less secure cases (e.g. B3 1340+439 and 
B3 2349+410 which are detected at two wavelengths only). 
The success is due to the fact that the fluctuations in m\ 
help in constraining some of the parameters, such as the 
difference between the two RMs and the intrinsic frac- 
tional polarisation mo. For instance, undetected points 
are taken as deep minima in m. The level of success is 
definitely lower (< 30%) in the second group. 

The results of the model have to be taken with cau- 
tion because of the limited number of data points. When 
reasonably secure, we give the average RM, the common 
CRM.obs and f c and the intrinsic fractional polarisation mo 
(vectorial sum of m(A = 0) of the two components). 

Comments on individual fits are contained in the notes 
to individual sources (|A.1|) . 



corrected to zero wavelength for RMq G . It has to be re- 
alized that RMq 6 has an uncertainty > 50 rad m~ 2 . 



